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Abstract 

We establish a characterization of coagulation-fragmentation processes, such that the induced 
birth and death processes depicting the total number of groups at time t > are time ho- 
mogeneous. Based on this, we provide a characterization of mean-field Gibbs coagulation- 
fragmentation models, which extends the one derived by Hendriks et al. As a by- product of 
our results, the class of solvable models is widened and a question posed by N. Berestycki and 
Pitman is answered, under restriction to mean-field models. 
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1 Introduction, objective and the context 

The time dynamics of a time homogeneous Markov process X(t), t > on a space Q = {77} 
of states r\ is described by the set of transition probabilities 

p^ V ;t):=F(X(t)= V \X(0) = C), C,V^Q, t > 0. 

Given the rates of the infinitesimal state transitions, the explicit expressions for the transition 
probabilities as solutions of a Kolmogorov system, are known only for a few special cases of 
the rates. The corresponding models are called solvable. For the above reason, time dynamics 
of Markov processes remain, generally speaking, a mystery. As an example, even for birth- 
death processes on the set of integers, the explicit solutions have been derived only for a few 
combinations of birth and death rates. This explains why the direction of research in this 
area turned to the estimation of the rate of convergence (=spectral gap) of the transition 
probabilities as t — > oo. Nevertheless, hunting for solvable models continues to be of interest. 

In the present paper we pursue the above objective for stochastic processes of coagulation 
and fragmentation (CFP's). We adopt the formulation of a CFP = CFP(N) given in [5] on 
the basis of classic works of Whittle [25] and Kelly [H] devoted to deterministic and stochastic 
models of clustering in polymerization, electrical networks and in a variety of other fields. A 
CFP Xff(t), t > is defined as a time homogeneous Markov chain on the state space Qn of 
all partitions n = (ni, . . . ,njy) : YliLi = N, of a given integer N. Here N codes the total 
population of indistinguishable particles partitioned into groups (=clusters) of different sizes, 
while ^ is the number of groups of size i. Infinitesimal (in time) events are a coagulation of two 
groups into one and a fragmentation of one group into two groups, and the basic assumption 
is that the rates (intensities) of the above two single transitions depend only on sizes of groups 
(and do not depend on N). Namely, the rate of a single coagulation of two groups of sizes % 
and j, such that 2 < % + j < N, into one group of size % + j is whereas the rate of a 

single fragmentation of a group with size i + j into two groups of sizes i and j is <p(i,j). The 
functions ip and <fi are assumed to be non negative and symmetric in 

Next, we define the induced rates of infinitesimal state transitions. Given a state rj e 
with rij, rij > for some 1 < i, j < N, denote by ij^'fi G f2jv the state that is obtained from r\ 
by a coagulation of any two groups of sizes % and j, and denote by K(rj — > rf 1 '^') the rate of 
the infinitesimal state transition n — > rj^^'. Similarly, for a given state r\ € f2jv with ni + j > 0, 
let T7(jj) be the state that is obtained from 77 by a fragmentation of a group of size % + j > 2 
into two groups of sizes i and j, and let F(r] — > T]nj\) be the rate of the infinitesimal state 
transition 77 — > r)(i,j). We assume that the rate K(rj — > 7] {i ' j) ) is equal to the sum of rates of all 
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single coagulations of n« groups with size % with rij groups with size j, and that F(rj — > f](ij)) 
is the sum of rates of all single fragmentations of n i+ j groups with size i + j into two groups 
of sizes i and j. As a result, we get the following expressions for the rates of state transitions: 

K( v ^rj^) = runrfiij), i^j, 2<i+j<N, 
K(n ^ riM) = Ui{Ui ~ ^ 2<2i<N, 

F (v^V(i,j)) = n i+j tj)(i,j), 2<i + j<N. (1.1) 

We note that an interpretation of the coagulation kernel K in terms of the kinetics of droplets 
of different masses can be found in [T9] . 

Following [IT] , we call CFP's with rates of state transitions of the form ( 11.11) mean field 
models, meaning that at any state rj £ Q N , any group can coagulate with any other one or can 
be fragmented into any two parts. We also note that ([5]) a characterization of positive rates 
of single transitions ip(i,j),<fi(i,j) that provide reversibility of mean- field CFP's is known 
(0). 

We now describe briefly the context of the present paper. The paper is devoted to the 
time evolution of the above mean field CFP's and it consists of two sections. Section 2 is 
divided into three subsections. In Subsection 2.1 we characterize the CFP's Xjsr(t), t > 
having time homogeneous processes |Xjy(t)|, t > depicting the total number of groups at 
time t > 0. The key result of the paper, stated precisely in Theorem 1 in Subsection 2.2, 
establishes the equivalence of the following two conditions: 

(i) The birth and death process |Xn-(£)|, t > is time homogeneous; 

(ii) The conditional distribution of a CFP(N), given a total number of groups at time 
t > 0, is a time independent Gibbs distribution on the set of partitions of iV with given number 
of components. 

Consequently, a characterization of Gibbs CFP's, which extends the one by Hendriks et 
al ( [13]), is derived. 

In Subsection 2.3 we discuss the following three topics related to our main result: Steady 
state distributions of CFP's, Gibbs CFP's on set partitions and Spectral gaps of Gibbs 
CFP's. In particular, under restriction to mean-field models, we obtain a negative answer 
to a question posed by N. Berestycki and Pitman ([3]) about the existence of certain Gibbs 
CFP's. 
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2 Main result 

We say that states fj,ij G are neighbors: fj ~ rj, if one of the states is obtained either 
by a single coagulation or a single fragmentation of components of the other state. Then the 
preceding description of a CFP(N), say (t), t > 0, starting from an initial probability 
distribution p on Vt^ allows us to write the corresponding Kolmogorov system as follows 

^p p (rj',t)(K(r}^ri) + F(Tj^ri)), f}, V e^ N , t>0. (2.2) 

Note that the seminal system of Smoluchowski equations (1918) for pure coagulation can 
be viewed as an approximation to (12. 2p obtained by neglecting correlations between group 
numbers at time t > 0. This issue is widely discussed in the literature, (see [5], [2D] , [8],[I]). 

groups. 

{riiit), . . . , n N {t)) G VL N) t > 0, a central 

TV 

ixjftt)! :=£>(*), ( 2 - 3 ) 
i=i 

which depicts the total number of groups in the generic CFP at time t > 0. We denote 
throughout the paper 

^N,r = {v e Q N ■ \v\ = r }> r = 1, . . . , iV 

the set of all partitions of TV with exactly r components. 

It follows from the definition of a CFP(N) that * > is a Markov birth and 

death process on the state space {1,2,..., iV}, with rates of birth and death 
K,n, 1 < r < TV — 1 . fi r ,N, 2 < r < AT, respectively, defined in a usual way, as in (I2.4p below. 
However, in contrast to the generic CFP(N), the process t > is, in general, not 

homogeneous in time, which presents a big problem for the investigation of the process. 
The following example demonstrates the phenomenon of dependence of the rates \ r ,N,Hr,N 
on time t > and on an initial distribution p of the generic CFP(N), that causes the time- 
inhomogeneity of the induced birth and death process. 

Example: Consider a CFP(N), N > 4 of pure coagulation with ?/>(l, 1) = ■!/>(!, 2) = and 



2.1 Process of the total number of 

In our study of time dynamics of a CFP X^\t) = 
role is played by the induced stochastic process 



transient 



6 



all other ip(i,j) > 0. It is clear that f2jv iA r_ 2 = 771 {Jrj2, where rji = (N — 3, 0, 1, 0, . . . , 0), 
V2 = {N — 4, 2, 0, . . . , 0). Assuming that the process starts from an initial distribution p on 
&n,n-2, s.t. p(rji) =pi>0, i = 1, 2, p 1 + p 2 = 1 and denoting A = J2ten N , N _ 3 K (Vi > 
0, % = 1, 2, we have 

PpiVi'it) = -AiPpiVht), t>®, 2=1,2, 
since the transitions f2/v,jv-i — > 77^, i — 1,2 are impossible. Hence, 

P P (Vi;t) = Pl e~ At , t>0, z = l,2 

and consequently, it follows from (12. 8p below that the rate of death 

A ie - Al Vi + A 2 e- A2 V 2 
m-2,N{t'i P) = zjTt 1 _ A , t , < > 

depends on p and £ iff A x = (N - 3)^(1, 3) ^ A 2 = ip(2, 2). ■ 

We further assume that the CFP's considered are ergodic. 

We first distinguish CFP's(N) which induce time homogeneous processes 

t > 0, under any initial distribution p on fi^r, i.e. processes with birth and death 
rates not depending on t > and p. Let A r) jv(t; p), 1 < r < N — 1 and p, r> N(t] p), 2 < r < N 
be, respectively, the rates of birth and death at time t > 0, of some : 



K,N{t', p) = lim 



p(|xW(* + At)|=r + ll 



At^0+ At 



p(|JfW(t + At)|=r-ll 



A»r,2v(*; p) = lim — . (2.4) 

( 12.41) tells us that under the assumption of ergodicity of the generic CFP(N), the time inde- 
pendence of birth and death rates implies their independence on p. 

Clearly, the birth and death rates in (12 .4p are implied respectively, by the rates ifi of 
single fragmentations and the rates of single coagulations of the generic CFP(N). It turns 
out that the required necessary and sufficient condition of time homogeneity of the process 
t > has a simple probabilistic meaning. 

Lemma 1 \X$(t)\, t > is a time homogeneous birth and death process under any initial 
distribution p on Q^, if and only if the generic CFP(N) is such that for a given 1 < r < N 
the sums of rates ^2fj^ v K(r] — > 77), r\ 6 &>N,r an d X^^-^X 7 ? — * V)> V e ^N,r do not depend on 
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rj G fijv,r- Under the above condition, the first sum and the second sum are equal to the rate 
of death p r ^ and to the rate of birth X r) ^ respectively, so that for any r\ e f£jv>> 

X r , N = hm + ^ P (\X%\t + At)\=r + l\ X%\t)=rj), l<r<N-l, 

Vr, N = Km + ^F(\X < g\t + At)\ = r-l\ X^\t) = r^, 2<r<N, (2.5) 

under any initial distribution p on Qn and all t > 0. 

Proof Recalling that t > is Markov, by the Markovian property of the generic 

CFP(N), we firstly assume that |X^(t)|, t > is time homogeneous, so that for all p on 
Qn and all t > 0, \ r ,N(t; p) = X r ,N, 1 < r < N — 1 and p r ,N(t; p) = p r ,N, 2 < r < N. We now 
rewrite (12. 4p as 

1 J2 v eu N P(\^\t + At)\=r + l I X^(t)= V )p(x^\t)= V 
\ r , N = hm — — — 

Kr < N- 1, 



// r jv = lim 



p(|X«(t)|=r 

1 E^P^Jftt + At^r-l I ^ ) (t)=r ? )p(x^(0=^ 



p(|JfW(t)|=r 



A^0+ At |> / I v 1 '"/ 

2 < r < iV. 

(2.6) 

By the ergodicity and the time homogeneity properties of the generic CFP(N), the limits 
f b ( V ;r,N):= lim i- P f \X%\t + At)\ = r + 1 I = A l<r<JV-l, 

/ d (r/;r,iV):= A lim + ^p(|X^ ) (t + At)|=r-l I X%\t)=ity, 2<r<N, (2.7) 

do not depend on t > and p, for all n 6 fijv,D so that 

/ b (r/;r,iV)p(x^ ) (t)=r ? ) 

Ar,iv= V ^— f r S i<r<7V-l, 

v#£,r P(\X^(t)\ = r) 

f d ( V ;r,N)p(x^\t)= V ) 

Hr tN = V ^— f r S 2<r<iV. (2.8) 

P(\X%\t)\=r) 

Next, setting in (j23j) . t = and p(C) = 1, for a C G ^n,v, (so that P^X^(O) = <f) = 1) it 
is easy to conclude that (12. 8p together with the time homogeneity assumption imply 

K,n = /b(C; r , = const, Pr,N = fd(C, r, N) = const, 
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for all C G f2jv,n which proves the necessity of the condition (12.51) . The sufficiency of 
(12. 5p follows immediately from (12.81) . after we observe that in view of (12. 7p . the quanti- 
ties fd(r];r,N), f b (n;r,N) are equal respectively, to the sum of rates of single coagulations 
J2^ v K(n — > fj) and to the sum of rates of single fragmentations F(rj — > fj) at a state 

T) G &N,r- ■ 

In the rest of this subseqtion we will treat the case when the time homogeneity condition 
in Lemma 1 holds, writing simply |X/v(£)|, t > 0. Now our objective will be to characterize 
the rates ip(i,j),(j)(i,j) that provide the condition (12.51) . The condition (12.51) says that for 
given N and r each one of the two limits in the RHS of (I2.5P is the same for all n G fijv.r and 
all p on Qjy. Consequently, the above condition conforms to two separate systems of linear 
equations, one for ip(i,j) and one for (f)(i,j), and each one consisting of \flN,r\ equations for 
each 1 < r < N. It is easily seen that for a fixed N there is a variety of solutions to each of 
these systems, which are valid for all 1 < r < N. 

For example, employing the aforementioned meaning of the limits fj, and fd, one can verify 
that for a given N > 3 the following rates depending on iV satisfy (12. 5p : 

and 

where l\ and l 2 are arbitrary nonnegative functions. 

However, due to our basic assumption that the rates ip and <fi do not depend on N, time 
homogeneity of the process |Xjv(t)| t > implies a very special form of the above rates of 
singular transitions. 

Proposition 1 {|Xjy(t)|, t > 0}n>i is a sequence of time homogeneous birth and death 
processes induced by a sequence of {CFP(N)}n>i with rates of single transitions ip(i,j) and 
4>{h ])-, if an d on ly if the above rates are of the form: 

= a(i + j) + b, i, j > 1, a > 0, 2a + b > (2.11) 

and 

v k := Yl <f>(iJ) = <f>0-Mk-l), k>l, (2.12) 

l<i<j: i+j=k 



i+j, if 2 < % + j < N - 1 
h(N), if i+j = N 



(2.9) 



0, if 2 < i + j < N 

l 2 (N), if i + j = N, 



(2.10) 
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where v(k) is the sum of rates of all possible single fragmentations of a group of size k > 2 
into two groups, whereas v% = 0(1, 1) > is arbitrary. 

Proof: We employ the preceding lemma. Assuming that the processes |Xjv(t)|, t > are 
time homogeneous for all A > 1, we apply the second part of (12. 5p with r = 2 to obtain 

^{i, N — i) = imp, i = 1, . . . , N - 1, N > 1. 

Therefore, 

il>(i,j) = s(i + j), i,j>l, (2.13) 

where s is some nonnegative function on integers which are greater or equal to 2. 
Next, consider the two states 771 , 772 G &n,3, N > 5 : 

nn-2 

■m = (2,0,. .. ,0,^,0,. ..,0), 

njv-3 

r/ 2 = (1,1,0,...,0,^,0,...,0). (2.14) 

Applying the equation f d (r)i] 3, N) = f d (rj 2 ] 3, iV), gives 

2^(1, iV - 2) + ^(1, 1) = V(l, iV - 3) + ^(AT - 3, 2) + ^(1, 2), (2.15) 

which by virtue of (12.131) . is equivalent to 

2s(N - 1) + s(2) = s(N - 2) + s(N - 1) + s(3), N > 5. 

Taking into account that the last relation should hold for all iV > 5, we rewrite it as s(k) — 
s(k — 1) = s(3) — s(2), k > 3, which proves the necessity of (12. lip . For the proof of the 
necessity of (12.121) we consider the quantities fb{v\ 2, N) for fixed and all states n of the 
form 

i N-i 

?7 = (o, . . . , o, . . . , 0,^, 0, . . . , 0,^, 0, . . . , 0) g n Nt2 , 1 < % < n - 1. 

Using the notation in (12.121) . the condition that fb(rj;2,N) should be the same for all the 
above r\ can be written as 

v(i) + v(N -i) = const, 1 < i < N - 1, (2.16) 

or, equivalently, v(N — 1) — v(N — 2) = u(2) — f (1) = u(2). Since the latter relationship should 
hold for all iV > 2, it implies ( I2.12p . We turn now to the proof of sufficiency of the conditions 
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( 12. lip and (12. 12[) . Supposing that (12. lip holds, we have for a state r\ G Qn,t '■ 

Uiini - 1) 



l<i<j<N l<i<N 



[1,1 Hi 



^(^>i n i - E ^ 

l<i,j<N l<i<N 



l<i,j<N l<i<N 



= - (2aNr + fer 2 - 2aiV - 6r) , r = 2, . . . ,N, 
f b (r);r,N) = ^ v ( k ) n k 

l<k<N 

= Y v(2)(k-l)n k = v(2)(N -r), r = 1, . . . , N - 1. (2.17) 



Kk<N 



Corollary 1 The rates of death and birth of a time homogeneous Markov process 
\XN(t)\, t > are given by 

Li r ,N = Z ( 2aN + rb\ 2 < r < N 



2 V 

A ri 7v = 0(l,l)(iV-r), l<r<N-l. (2.18) 

Remark 1 

(i) The birth and death process \Xj^(t)\, t > with rates given by A2.18\) . has the following 
interpretation, not related to the generic CFP(N). Consider a nearest neighbor spin system 
(for reference see fT7\l ) of "0"-s and "l"-s on a complete graph on N vertices (sites). Assume 
that one of the sites is occupied with a " 1" which never flips, while spins at all other sites 
perform flips — > 1 and 1 — >• with rates X r ^N and p, r -i,N respectively, where r is the total 
number of sites of the graph occupied by "Vs. (The latter says that a site occupied by a 
"1" has r — 1 neighbors occupied by Is and a site occupied by a "0" has r such neighbors). 
Consequently, at a state with r > 1 "V s, the total rate of — ► 1 flips is X r> N '■= (N — r)X r ^ 
and the total rate of 1 — ► flips is fi r> N '■= (r — l)p, r -i t N- Therefore, the induced birth and 
death process, say Civ(^), t>0, on {1, . . . ,N} depicting the number of sites occupied by "Vs 
at time t > is Markov and time homogeneous. Clearly, if 

K tN = 0(1,1), l<r<N-l, 



-foaN + brj), 2 < r < N, 
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the process Ov(t) conforms to the process |Xjv(£)|, £ > 0, associated with CFP's(N) given by 
li2.11\) . F2.12\) . Finally, it is appropriate to note that after interchanging the roles o/"0"-s and 
"V'-s, the spin system with the rates X r , p r -i as above, is known (for N fixed) as a contact 
process. 

(ii) It follows from Proposition^ that the class of CFP's(N) that induce time homogeneous 
processes |Xjv(t)|, t > includes processes of pure coagulation (<fi{l,l) = in 112.12]) ) and 
processes of pure fragmentation (a = b = in A2.11\) ). Also observe that the time homogeneity 
requirement determines uniquely the form of rates of singular coagulations, while it leaves a 
certain freedom in the choice of rates of singular fragmentations. 

As far as we know, there are no explicit solutions, i.e. explicit formulae for transition 
probabilities P(|Xjy(t)| — r), t > 0, 1 < r < N, for birth-death processes with the rates given 
by (12.181) . when a, b > 0, 1) > and the initial distribution is concentrated on some state 
C G The problem here is that the birth and death rates in (I2.18P are polynomials in r of 
different degrees, which are 1 and 2 respectively. A survey of solvable birth-death processes 
with polynomial rates is given in [22J. 

We will see in the next subsection that under the above condition of time homogeneity of 
I^jv(^)I) t >0 and certain initial distributions p, the corresponding CFP's(N) are solvable. 



2.2 Solvable CFP's 

Let a CFP(N) considered start from an initial distribution p on Qn, with projections p r on 
the sets flN,r, r — 1, . . . , N: 

I p(r] I \rj\ = r), 7iett N , r , if p(tt Ntr ) := p(\rj\ = r) > 
Pr(v) ■= < (2-19) 

I 0, 7] G Qjy,r , if P{^N,r) = 0. 

It is in order to note that the set of all distributions p on with given projections p r , is 

N 

{p : p(r]) = p r (rj)p(Q N>r ), r}eVL N) r, y^ y p(Q N>r ) = 1, p(Q N ,r) > 0, r = l,...,N}, 

r=l 

(2.20) 

i.e. the projections p r , r = 1, . . . , N define the associated distribution p up to the factors 
p(p N ,r), r = l,...,Nia d22Q]). 
Accordingly, we write 



P P (v,t) = w[xj>\t) = v I \x%\t)\ = \v\) v{\x<£\t)\ = \ v \ 

:= Q{ V ,p;t) b{\r]\,p;t), V eQ N , t>0, (2.21) 
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where Q(r], p; t) and b(\rj\, p; t) denote respectively the first and the second factors in the RHS 
of the first line, while the conditional probability Q obeys the initial conditions 

Q{v,p;0) = p r (v), V e n N>r , r = l,...,N, (2.22) 

for any p on fl^, which follow from the definitions of Q and p r . 

We will be interested in CFP's possessing a conditional probability Q not depending of 
time under certain initial distributions p. If this is the case, it follows from (I2.22p that 

Q(rj, p; t) =pr(v), V£ ^N,r, r=l,...,N, t>0. (2.23) 

Obviously, time independence fl 2 . 2 3 [) holds for any CFP(N) starting from its stationary dis- 
tribution. In view of this, we adopt the following convention. 

Definition 1 A CFP(N) possesses a time independent conditional probability Q if (12.231) 
holds for certain projections p T on f)jv r , t = 1, . . . ,N and all initial distributions p on fl^ 
from the associated set (12.201) . 

Next we write the Kolmogorov system for a birth and death process \X^ (t)\, t > 
with rates A rj jv(i; p), p r ,N(t] p) '■ 

b(r, p; t) = -b(r, p; t) (K,n{^ P) + Pr,N(t; p)) + b{r + 1, p; t)p r+ltN (t; p) + 

b(r-l,p;t)X r -i, N (t;p), r = l,...,N, (2.24) 

where 6(0, p; t) = b(N + 1, p; t) = 0, t>0. 
The following assertion is crucial for our study. 

Proposition 2 The following two conditions (i) and (ii) are equivalent. 

(i) A CFP(N) possesses a conditional probability Q independent of time t > 0; 

(ii) The birth and death process \Xr?(t), |, t > is time homogeneous. 

Moreover, the projections p r , r = 1, . . . , N defining by 112. 23\) the time independent condi- 
tional probability Q are the unique solution of the two systems of equations: 

Pr+l,Npr(v) = ^2 Pr+l(C)K(( ~^ V), V € ^N,r, r = 1, . . . , JV - 1, (2.25) 

K,NPr+i(0 = Pr(v)F(v^0, Ce n N , r+li r = l,...,N-l, (2.26) 

where the rates of state transitions F and K are given by M.l\) . 112.11}) , 112.12}) , while the rates 
of birth and death are as in 112.18}) . 
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Proof: We substitute (12.211) in the Kolmogorov system (12.21) to obtain 
b(r, p; t)Q(7], p; t) + Q(r], p; t)b(r, p; t) = 

-Q( V ,p;t)b(r,p;t)( ^ K( V ^() + E F ^ ~+ 0) + 

6(r + l,p;t) £ Q(C,p;*)ir(C->77) + &(r-l>P;*) ^ ^CtfWC - 
rj e Q N , r , r = 1, . . . , N, t>0, 

{l ,N = Vn,n + i = 0, b{0,p;t) = b{N+l,p;t) = 0, t > 0. (2.27) 

We firstly prove the implication (ii) =>- (i). We substitute in the LHS of (12.271) the expres- 
sion for the derivative b(r,p;t) from (12.241) . assuming that the process |Xjv(t)|, t > is time 
homogeneous and that the initial distribution is p. Then, by virtue of Lemma [TJ the system 
(gj27D becomes 

Q(r), p\ t) (b(r + 1, p; t)fj, r+ljN + b(r - 1, p; t)A r _ ljA rJ + Q{rj, p\ t)b(r, p; £) = 

6(r + l,p;t) Y, Q((,P;t)K((^r,) + b(r-l,p;t) ]T Q(C, p; t)F(( - 77), 

T]eQ N , r , r = l,...,N, t>0, 

Q N = fijv,jv+i = 0, 6(0, p; t) = fe(iV + 1, p; t) = 0, t> 0, for all p on f^jv, (2.28) 

where by the assumption made, the rates K(( — > rj) and F(rj — > £) of state transitions are 
implied by (12.111) . (I2.12p respectively and the birth and death rates are as in Corollary 1. 
Given an initial distribution p, constants a, b and fragmentation rates <f){i,j) obeying (I2.12p . 
a finite Kolmogorov system ("12. 28[) has a unique solution Q, provided a 2 + b 2 + 0(1, 1) > 0. In 
particular, (I2.28j) is satisfied by the time- independent Q, such that (12.231) holds for all initial 
distributions p with the projections p r , r = 1, . . . , N that obey (12. 25p . (12.261) . 

To prove the implication (i) =>■ (ii), we observe that by virtue of (12. 8p . the condition (i) 
implies that the birth and death rates do not depend on t > 0. By our remark after (12.41) this 
leads to the conclusion that the rates do not depend on p either. 

Next, we set t = in (I2.28P to derive by virtue of Definition 1 and the time homogeneity 
of |Xjv(t)|, t > 0, that the projections p r , r = 1, . . . , N should obey (12.251) and (12.261) . 

It is left to show the existence and uniqueness of the solution p r , r = 1,...,N for the 
system of equations (I2.25p . (12.261) . where the rates of state transitions are induced by tp and 
4> as in (12. lip . (12. 121) . Recalling Lemma [fl we treat the ratios 

KU -> rj) 

Pc{( -> V) ■= > C e &N,r+l, V e ^N,r, C ~ v, 

Pr+1,N 



transient 



14 



r = 1,...,7V-1 

as the one- step transition probabilities of a discrete time nearest-neighbor "coagulation" 
random walk on the set of partitions Qjy. Then p r (rj), n E ft N>r , m the set of equations (12.251) 
can be interpreted as the probability that the random walk starting at n* = (N, 0, . . . , 0) £ 
Qn,n reaches a given state r\ £ f)jv,r a t the (N — r) — th step, so that (* = (0, . . . , 1) is 
the absorbing state. In a similar manner, we consider the nearest neighbor "fragmentation" 
random walk on with the transition probabilities 

F(r] -> C) 

Pf(v -* := — > , v e Q N , r , C e fijv, r +i, C ~ »7, 

r = l,...,iV-l, 

that starts at (* = (0, ... ,0, 1) £ Qn,i- In this case, p r (rj), rj £ J7jy,r in the set of equations 
(I2.26P is the probability that the "fragmentation" random walk reaches a given state n £ Ojv, r 
at the (r — 1) — th step, r/* = (TV, . . . , 0) being the absorbing state. Clearly, each one of the 
two systems has a unique solution p r , r = 1, . . . , N whenever a 2 + b 2 > in the first case and 
0(1, 1) > in the second case. 

We demonstrate that when (a 2 + 6 2 )0(1,1) > (=both coagulation and fragmentation 
hold), the two systems of equations have the same solution if and only if the transition prob- 
abilities Pq and Pp are related in the following way. Let p r , r = 1, . . . , iV be the probabilities 
corresponding to the "coagulation" random walk, under some fixed a, b : a 2 + b 2 > in 
(12. lip . Then the equations (I2.26P for the "fragmentation" random walk have the same solu- 
tion p r , r — 1 , . . . , N if and only if 

p r {r])P F (n -»• C) = p r+ i(()P c (( -»• rj), T]E Q Ni r, C e fitf,r+i, V ~ C- (2-29) 

The sufficiency of (I2.29P is seen immediately, while the necessity can be derived from the 
following general reasoning, based on the observation that each one of the equations (12.251) 
and (I2.26P is time reversal of the other one. Let p r , r = 1, . . . , iV be the common solution 
of (12.251) and (I2.26p . Then from (I2.25p . applied for r = TV — 2, we conclude that under given 
Pn-i and pn-2, the values a, b in (12.111) are uniquely determined. Hence, p r , r = 1, . . . , N 
uniquely determine all probabilities Pc in (12.251) induced by (12. lip . If now some Pp satisfies 
(I2.26P under the above p r , r = 1, . . . , TV, then (I2.25P should be satisfied by 

Pr+l(() 

The aforementioned uniqueness of the probability Pq proves the claim. (In the discussion 
following the proof we find explicitly the solution p r , r = 1, . . . , N and demonstrate that the 
rates of singular fragmentations derived from (I2.29p satisfy the condition (12.121) ). ■ 
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Our next purpose will be to find explicitly the solution p r (r]), n e fijv,n r = 
of (I2.25p . (l2.26p . in the case when the time homogeneous process |Xjv(£)|, t > is given by 
(I2.18p . Since the sets £In,i,^n,n are singletons, it follows from the definition of the conditional 
probability p r that Pn{v*) = 1> Pi(C*) = 1- The following two cases should be broadly 
distinguished. 

Case 1: Non zero coagulation, i.e. a 2 + b 2 > 0. 
Following the illuminating idea of Hendriks et al ([13]), we will seek the probabilities p r in 
question in the form 



Pr{v) = PN,r(v) = [B N> rJ 



n\\ni\ ■ ■ ■ tin- 

•n=(n 1 ,...,n N )eQ Ntr , r = l,...,N, (2.30) 

where is the normalizing constant (= partition function) known as the (N, r) partial Bell 
polynomial (see e.g. [I], [21]) induced by the sequence of weights {a k }^° that do not depend 
neither on iV nor r. It follows from (12.301) that for given n = (rix, • • • , n^) G fijv, n such that 
n i+ j > for some 2 < % + j < N, and ( = J](ij) G flN,r+i, 



Pr+l(C) = ( B N ,r \ j (I'Z ) ( ) ' lf * ^ J 

p r {rj) \B Ntr+1 J } (±\ ( m \ ifi=j 

Hence, setting, in accordance with Proposition [1] and (11.11) . 



(2.31) 



(a(t + j) + b)(n t + l)( nj + l), xii^j 
K (V(i,j) -> V) = { . , n ., , (2-32) 

(2m + fe) (n ' +1) 2 (w ' +2) , otherwise, 

where a 2 + fr 2 > 0, the equations ( 12. 25ft conform to 

% \ A(aA: + 6)E 



Pr+1,N — I "5 : — j >^ ~ Kk, 

(ni, . . . ,njv) G fijv,D r = 1, . . . ,iV — 1. (2.33) 



Since the RHS of (12.331) should not depend on 77 e f2/v,r the equations are solved by the 
weights defined recursively by 

(ah + b) Y". , ■ , a,- a,- 

I ^l±^fe 1 fc>2 L (234) 

2(fc — 1) ~ V ; 
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This is just the solution obtained, by quite different considerations, in [13J (see (18) there), 
for pure coagulation processes. 
Continuing (12.331) . we get 

B N 

fJ>r+l,N = ) VVfc - l)w fc , (2.35) 

(ni, . . . ,n N ) G fij^r, r = 1,...,N — 1, (2.36) 

which leads to the following relation between the constants u. r+liN , B Njr , B Nr+ i induced by 
the weights (jZ33j) : 

// r+1 , w = (iV-r)f-^-), r=l,...,tf-l. (2.37) 



N 

Taking into account that Sjv,jv = = (^0 _1 > we S e ^ the explicit expressions for the Bell 
polynomials in the case considered: 

B ^ = ^0^V r = l,...,JV-l, (2.38) 

where fJLi^ as in (12.181) . Remarkably, the expression ( 12. 38ft for the Bell polynomials enables 
us to find explicitly the weights a^, k > 1, without solving the recurrence relation ( 12. 34ft . In 
fact, by (12T30D . 

T-riV 

which can be written as 

T\ k (k a +—) 

ai = 1, a fc = llr - 2 V; — , *: > 2, 2a + 6 > 0, a > 0. (2.40) 



Remark 2 TTie recurrence relation {2.3$ can be viewed as a modification of the classic con- 
volution formula, 



| ^2 a i a v k = 2,..., a 1 = l, 



i+j=k 



which determines the Catalan numbers (see e.g. fTSf ). It is interesting to find the generating 
function g(x) = YlT=i a k xh f or ^ e sequence of weights defined by {2.34 )• Setting 



y(x) = it follows from {2.3$ that the function y obeys the differential equat 

n i 2a + b 

y (1 - axy) = y a 2 , a 2 = — - — , y{0) = at = 1, 

which implicit solution is given by 

7 „ 2a + b 



ton 



y(x) = (l + -xy^j b , b > 0. 
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We now recover the fragmentation rates given by f!2.29j) in the case of coagulation rates 
fl2T32|) . By fT2~311 . fl2~37l) and f[2~T8D we have 

h ■' ) ( a(i +j) + b), if % ± j 
F( V ^ = 0(1,1) { I ^ S J} r rJ (2 .41) 



(2ai + bj, if i = j. 



for r] = (m, . . . , un) G flN,r, such that n i+ j > for some 2 < i + j < N, and £ = G 
f2jv,r+i- Note that by (12.341) . the rates of singular fragmentations induced by (12.411) satisfy the 
condition (12.121) . This latter condition appears to have a physical meaning in the context of 
CFP's describing polymerization (see [23J). 

Also note that in the case considered, p r (v) > 0, V £ ^JV,r> r = 1, ■ ■ ■ , N and that the 
rates of single transitions are determined by (12.411) up to the constants a, b and 0(1, 1). 

Case 2: Pure fragmentation. 
It is clear from the preceding discussion that under the fragmentation rates of the form (12.411) 
(with a, b that are not related to coagulation rates) the solution p r of (12.261) is given by (I2.30p . 
However, in contrast to the case of pure coagulation, Proposition [1] leaves freedom for the 
choice of rates of single fragmentations obeying (12.121) . In view of this, the probabilities p r 
solving (12.261) will depend on a particular choice of the above rates, so that p r , will be of 
Hendriks et al form (I2.30P if and only if the rates of single fragmentations are induced by 
(I2.4ip . This is illustrated by the toy example below. (We recall that under all above choices 
of rates of single fragmentations, the rates of the induced pure birth process remain the same: 
A r = </>(l,l)(N-r), r = l,...,N). 

Example Let 

{(i + j — 1), if % — 1 or j — 1 
(2.42) 
0, otherwise, 

so that (12.121) holds. The corresponding fragmentation random walk is in effect a deterministic 
chain on TV states (i, . . . ,(r, ■ ■ ■ ,(n, such that 

N-r+l 

( r = (r - 1, 0, . . . , 0, , 0, . . . , 0) G Sl Ntr , l<r<N-l, 

Gv=(iV,0,...,0). 

Consequently, (I2.26P implies p r iji) = ^-c,X r i)i V ^N,r, which is not of the form (I2.30p . 
The preceding discussion is summarized in our main theorem. 
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Theorem 1 Solvable CFP's 

Mean -field CFP's X^\t), t>0 with rates of single transitions $2.11\) and \2. 7j| ) and initial 
distributions p on Qn, s.t. p r on &N,r, r = 1, . . . , N satisfy Ii2.25\) . f2.26\) . have time dynamics 
given by 

Ppiv, t) = Pr{v) h i r -, p;t), v e r = i,...,N, t>o, (2.43) 

where b(r, p; t) are transition probabilities of the associated time homogeneous birth and death 
process \X^\t)\, t > with rates $2.18\) . In partucular, if rates of singular coagulations are 
positive, then p r is given by 112. 30\) . F2. 38\) and 1^2. 40\ ), while in the case of pure fragmentation 
p r , r = 1, . . . , N satisfying 112.26]) are as before, if and only if 1^2. 4 1\) holds. 

Note that under b = in (|2.18p . the corresponding birth and death process is known as 
the Ehrenfest process (=urn model). 



Remark 3 .Initial distributions p. 

CFP's(N) with single transitions Ii2.11\) . f2.12\) but with initial distributions p that do not 
obey the equations A2.25\) . 112.26)) in Proposition 2 are not solvable, since in this case the 



conditional probability Q depends on t > and p, though the processes \X^ (t)\, t > 
time homogeneous. 



are 



Remark 4 .Transition rule for Gibbs fragmentation. 

We now explain that the probabilities Pf induced by the positive fragmentation rates j2.41\) 



define the following simple rule of a state transition via a fragmentation from n £ &N,r to 
V(ij) e Sl N ,r+v By [273$ , 

p r( ,^ (ij . ))= (' + i-l>"^ H Tr)' 1 - ***** (2.44) 



2 Tl l l+m=2i a l a m) , tfi—j- 



N -r 

Under a given rj = (n\, . . . ,un) £ &N,r, the first factor in the above expression is the 
probability that a component of size i+j > 2 is selected to fragmentate, while the second factor 
specifies the probability that, conditioned on the first event, the selected component splits into 
two components of given sizes % and j . As a result, \2.4$ conforms to a transition procedure 
postulated in j^j, in which the first and the second factors are called the linear selection rule 
and the Gibbs splitting rule respectively. Theorem 1 and Proposition 2 say that the mean-field 
transition mechanism \2.44 ) is forced by the requirement that the process |JGv(t)|, t > is 



time homogeneous and the rates of singular coagulations are positive. 
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A historical note This note concerns exclusively the research on solvable CFP's. Time 
evolution of the stochastic model of pure coagulation was formulated by Marcus ([19"]) who 
was also apparently the first to reveal the relationship between Kolmogorov equations and 
its deterministic analog presented by Smoluchowski equations. Solutions to Smoluchowski 
equations for pure coagulation with kernels K induced by (f>(i,j) = const, <p(i,j) — i+j and 
<fi(i,j) = ij were obtained long ago by researchers in the field of colloid aerosol chemistry (for 
references see [T9],[I]). Lushnikov ([TS] derived explicit formulae for the expected numbers 
Enj(t), t > 0, j — 1, . . . , TV for the process X^\t), t > of pure coagulation with (f)(i,j) = 
i + j, hj > 1) with the help of the generating function for transition probabilities p^(n; t), t > 
0, rj G Ojv- The aforementioned stochastic model is known as the Marcus-Lushnikov process. 
In ([IB], treating Smoluchowski equations as an approximation to Kolmogorov ones, Luchnikov 
proved the important fact that the solution to Smoluchowski coagulation equations with a 
general coagulation kernel, can be presented as a mixture of Poisson distributions with time 
dependent parameters. (Note that these parameters were found explicitly for the Marcus 
-Lushnikov model only). A further important contribution was made by Hendriks, Spouge, 
Eibl and Schreckenberg ([13]) who found explicitly the transition probabilities Pc_{n\t), t > 
0, C,n e Qn for a more general Marcus -Lushnikov model with (j)(i,j) as in (12. lip . This 
result, proven via a combinatorial argument, is based on the representation (I2.2ip with time 
independent conditional probability Q. 

2.3 Discussion of the main result 

• Steady state distribution. Firstly, consider solvable CFP's with nonzero rates of 
single coagulations and fragmentations. By (II. ip . (I2.32p . (12.411) and Theorem 1, the 
implied rates i/j, 4> of single coagulations and fragmentations respectively, are 





where a-,-, j > 1 are given by (12.401) . Thus, the ratio of the above rates is equal to 



ij(i,j) _ aittj 



i,j>l. 



(2.45) 



(f>(i,j) 0(1, 



■i+j 



Setting in (I2^5jl a* 



#1,1) 



gives 



if>{i,j) _ cucij 



i,3>l 



<f)(i,j) a i+j 
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which shows that the criteria of reversibility of mean- field CFP's (see [5j) is fulfilled. 
Moreover, by Theorem 1, the above process is the only reversible process, within the 
class of solvable mean-field CFP's. By virtue of (I2.43p . the invariant measure of the 
process considered is 

u N (ri) = b{r) Pr {rj), rj G Q N>r , r = 1, . . . , N, (2.46) 

where b(r) = lim^oo b(r, p; t) is the invariant measure of the associated ergodic birth 
and death process. The probability measures p r on VLn,t, r = l,---,N defined by 
f )2.30p . (12.381) . (12.401) . belong to the class of multiplicative measures (= Gibbs distribu- 
tions) which play also a role in the theory of random combinatorial structures (see 
[21], [21], [6]- [8]). The explicit expression for the measure un is obtained in a straight- 
forward way from the known form (see [5]) of the invariant measure of a reversible CFP 
with rates fl2T4"oD : 

~rti~ri2 7, n N 
-1 ( a l a 2 ■ ■ ■ a N 



vn(v) = M , ,"• r,eQ N , (2.47) 

where cn is the partition function of the measure v^. Next, we embark on an analysis 
of the asymptotic behaviour of the measure u^, as N — ► oo. For this purpose we need 
to know the asymptotics of the weights {ak}^- By (12 .40 p . 

' b \fe-l 



^nt 2 (f^+o = (i) fc - 1 (^(^)(^ + i))~ 1 (^), , if mo 

if b = 0, 



! k\ u - z ■ n ■ it ■ / \ n 1 . ■ .... , .... 

a k =l \ /k+i (2.48) 



k — 1 fcfe — 1 

~~kT 



where k > 1 and (z) n := z{z + 1) . . . (z + n — 1) = T ^J^ is the Pochhammer symbol. 
Applying the Stirling's approximation, gives, as k — ► 00, 



dC^k-2, ifa6>0 



ak ~ < 



-'2 



if a = 0, 6 > (2.49) 



3 



C 3 C| jfe-f, if 6 = 0, a > 0, 

where Ci = Ci(a, 6), C 2 = C 2 (a, 6), C 3 = C 3 (a), C 4 = C 4 (a) are positive constants. The 
measure in (I2.47P is invariant under the transformations of the weights a k — > C k dk, 
with any constant C > 0. Thus, the asymptotic behaviour of the measure un considered 
is identical to the one with the weights 

d Art, ifa6>0 



a' h ~ < 



const, if a = 0, 6 > (2.50) 
C 3 Art, if 6 = 0, a > 0, 
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as k — > oo. In accordance with the classification suggested in [3] for multiplicative 
measures with regularly varying weights ~ k a , k — > oo, the measure considered 
belongs to the convergent class (a < — 1) in the first and the third cases in (12.501) . while 
in the second case in (12.501) it belongs to the expansive class (a > — 1). It was shown in 
[3] that the convergent class of z/jy exhibits a strong gelation, as N — > oo: with a positive 
probability all groups cluster in one huge component of size close to N. In contrast to this, 
(see [8]), the expansive measures v N have, with probability 1, as N — ► oo, a threshold 
value iV«+2 for the size of the largest group in the associated random partition. In the 
context of the CFP considered the above crucial difference is easily explained by noting 
that the first and the third cases in (I2.50p correspond to a "strong" coagulation, while 
the second case corresponds to a coagulation with a constant rate. 

Clearly, pure coagulation and pure fragmentation processes Xjsf(t), t > have the 
absorbing states rj* — (0, . . . , 1) and C* = (N, 0, . . . , 0) respectively. 

In the conclusion consider a non- solvable CFP(N) as in Remark 3. In view of the 
ergodicity of this process, its invariant measure will be identical to the one of a solvable 
CFP, starting from any distribution p on fi/v with the gibbsian projections p r on fijv r , 
r=l,...,N, given by (T2T30|) . 

• CFP's on set partitions. 

These are processes with values in the space Q[n] of partitions of the set [N] = {1,2,..., iV} 
(=set partitions). From the physical point of view, this means that, in the setting of 
this paper, the N particles are labelled, so that clusters forming a state of the process 
are subsets of the set [N]. Q\m- valued CFP's are a generalization of Kingman's co- 
alescent that provided a mathematical framework for a variety of genetic models, in 
particular the Ewens sampling formula. Kingman's theory, which is surveyed in |21j, 
is based on the theory of exchangeable partitions. The development of Kingman's co- 
alescent by Pitman [21] and his colleagues lead to Gibbs partitions as distributions of 
fl\N]~ valued irreversible processes of pure fragmentation or pure coagulation. Formally, 
the linkage between Q^- valued and Q[n]- valued CFP's is expressed via a simple com- 
binatorial formula and it is discussed in (21], [I], [7]J. Among CFP's on Q[n], Gibbs 
fragmentation processes introduced in [1] by N. Berestycki and Pitman play a central 
role. These processes are defined as time homogeneous Markov chains Il(t) e firjvi, t > 
of pure fragmentation, such that the conditional distribution of II (t) given the number 
of blocks of the random set partition Tl(t) is a microcanonical Gibbs distribution not 
depending of t > 0. In terms of CFP's on Q^, the above conditional distributions are 
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just the distributions (12.301) on flN,r, 1 < r < N. Correspondingly, the time reversal 
of the above process is called Gibbs coagulation. In [I] the authors posed a problem 
of characterization of the weights (in their notation) Uk '■= dkkl for which there exist 
Gibbs fragmentation processes, and they proved that, under the assumption that the 
fragmentation rates are defined by recursive and selection rules (12.441) . the unique Gibbs 
distribution is given by the weights (12.341) . In [I], p. 393 it was conjectured that some 
other, more complicated splitting rules might be of interest. We will demonstrate (see 
Proposition [3] below) that the aforementioned characterization of weights is valid for a 
broad class of fragmentation rules, that includes the above one in [I]. 
The problem reduces (see Problem 2 in 0]) to the characterization of weights a&, k > 1, 
(not depending of N) and transition probabilities of fragmentations Pp that satisfy 

: 

Ptv,,+i(C) = PnAv)Pf(v^0, (2-51) 

C e fijV,r+l, T = 1, . . . , N- 1, 

when p7v,r is a Gibbs measure (I2.30P on £In,v Regarding the probabilities Pp, we assume 
that they are of the following general form implied by the mean- field property: 

Pf(v - = ^f^'' 3 ^ V=(n 1 ,...,n N )e Q N , (2.52) 

where is a symmetric nonnegative function not depending of iV and 

c (v) = Si<i<i<AT n i+i^(*)i) ^ s ^ ne normalizing constant. Clearly, (12.441) is a particular 

case of fl232|) . 

For our subsequent considerations it is important to note that in (I2.30P all weights 
afc, k > 1 should be positive, due to the fact that = 1, N > 1. 

Proposition 3 Under the assumption \2. 5S\) . Gibbs distributions Pn,t, t = 1,...,N 
satisfy the equation $2.51\) if and only if the weights in $2. 30\) are given by \2. 31$ and 
the rates 4>(i,j) of single fragmentations are the same as in \2.J^1^ . 



Proof We assume that Gibbs distributions pN,r, 1 < r < N satisfy (12.511) . Treating 
(12.51)) when r = 1 and 

i N-i 

( = (o,...,^,...,o,^,o,...,o) en N , 2 , i = i,...,n-i, 



gives 



transient 



23 



b n,2\{ a N \ /<j)(i,N -i) 



ii N ^2i and 

B 2 i,2 \ ( a 2 i \ f<fi(i,i 



1, (2.53) 



B 2 i,i^lafJ\ v 2i ) 

ii N = 2i, where in both cases v k > is defined as in (I2.12p . Since Bn,i = a at, we have 



a i a N-i, if N 7^ 2i 



<(j) (i,N-i) = -^{ (2.54) 

¥1 



B m \ \al HN = 2i. 
Secondly, applying (12.511) for r = 2 with 

C e n N , 3 ■■ ((fa) = ((k 2 ) = ((N -fa- k 2 ) = l, 

where fa, k 2 , N — fa — k 2 are distinct positive integers, gives 

1 = E^4tT^-C), (2.55) 

where 

?7i e ^at,2 : 77i( fc O = Vi( N ~ fa) = 1, ??2 e ty/\r,2 : 772(^2) = ^(-/V - fc 2 ) = 1, 

% e ^,2 : 773(^1 + fa) = Vs( N ~ fa - k 2 ) = 1 

denote the three states from which it is possible to arrive, via one step fragmentation, 
at the above state ( G ^jv,3- Substituting in (12.551) the expression (I2.54p for and using 
(I2.3ip . we obtain 



Bn,3 \ f a N-ki v N-ki a 'N-k2 V N-k2 _|_ 

B N , 2 ' ^BN-k!,2(v kl + VN-ki) B N -k 2 ,2(Vk 2 + VN-k 2 ) 
a 'ki+k 2 v k 1 +k 2 



Bk 1 +k 2 ,2{Vk 1 +k 2 + v N-k!-k 2 ) 

We set now for a given iV > 3, 



(2.56) 



Mk) := . ^ r, 2<fc<JV-l. 

B k)2 {v k + v N -k) 

This allows us to rewrite (I2.56P as 

f N (N - fa) + f N (N - k 2 ) + f N (fa + fa) = C(N), N > 3, (2.57) 
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where C = C(N) is a constant w.r.t. fa, k 2 : N - fa > 2, N — k 2 > 2, fa + k 2 > 2. 
The solution of (12.571) is given by a linear function 

/jv(ife) = A^/c + B N >0, 2<k<N>3 

and the constant C = 2A N N + 3B N , N > 3, where the reals A N , B N : Ajy > 0, 2A N + 
Bn > 0. As a result, the following relation is derived 

- . akVk r = A N k + B N , 2<k<N>3. (2.58) 

B k ,2{Vk + v N - k ) 

We will show that (I2.58P forces the weights a k , k > 2 to satisfy (12.341) . Let 

< H k := lim sup (A N k + B N ), 

N~*oo 

for any fixed k > 2. H k = oo is impossible because v k > 0, k > 2, by (I2.54p . Hence, 
H k > is finite for all k > 2, which implies 

< A := lim sup A N < oo, B :— lim sup B N < oo. (2.59) 

N^oo N->oo 

Recalling that v x = 0, we apply f[23B]) with N = k + 1, k>2 and N = 2k, k > 2, to 
get 

= A k+1 k + B k+ i, k>2 

-Dfc,2 

and 

"'" -- A 2k k + B 2k , k > 2 



2-Bfc,2 

respectively. In view of (12.591) . the last two relations are in agreement if and only if 
A = B = 0, so that from (12.581) we recognize that 



Consequently, letting 

and denoting 

one obtains from (12.581) 



lim VN- k = oo, k > 1. 



V VN ^ 1 

z := lim sup > 1, 



a k v k 

= e k > 0, k > 2, 



B. 



k,2 



e k = lim VN- k (A]yk + B^) = z (ak + 6), k > 1, 

N^oo 
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where 



< a := lim v^A N < oo, b =: lim v^B N < oo. (2.60) 

N^oo iV-+oo 



Substituting the expression for e k into (I2.58P leads to the following relation 

z~ k (ak + b) 



v k + v N - k 

which implies 



A N k + B N , 1 < k < N - 1, 



^(a* + S) + ^»-»>(a<*-*) + 8) i<* <JV -i. (2.6i) 

Supposing z > 1, we obtain, in view of (12.601) . 

z~ fc ( s & + 6) = lim v N - k (NA N + 2Sat) = 
oo, if a ^ 

ife > 1. (2.62) 



z 



~ k ( \im.N-*oo(NvNAN) + 2bj , otherwise. 



In both cases this leads to contradiction, since in the case a = 0, we should have b > 0, 
by the definition of e k . Hence, z — 1. By CI2 .6 If) . this means that for a fixed iV, the sum 
Vk + VN-k does not depend on k, so that v k is linear in fc, namely v k = 0(1, l)(fc — 1), 
since i?i = 0, t> 2 = 0(1, 1). As a result, (I2.58P becomes 

afc( *r 1} = (iV - 2)Ajvfc + (iV - 2)5* = |*±* :=ak + b, 2<k<N>3, 

Bk,2 0(1,1) 

since the LHS does not depend on N. Recalling now the definition of B k 2 , gives (I2.34p . 



Remark 5 (i) In OV, Berestycki and Pitman characterized the Gibbs solutions of A2.51\) 
in the particular case of state transitions \2.44 )■ Our solution §2^JW derived under a 



more restricted mean-field assumption Ii2.52\) has the same form as in OV, but with less 
freedom on the constants a, b. 

(ii) The weights w k = a k k\ in the form of a finite product of linear factors appear also as a 
solution of a quite different characterization problem. Gnedin and Pitman JM/, extending 
Kerov's result fTh^ . proved that an infinite sequence {11*}, N > 1 of Gibbs random 
partitions of [N] is exchangeable if and only if in I12.3(J\) the weights, say w k = a k k\, are 
of the form 

fe-i 

w k = Yl (bl ~ aj , k > 2, w x = 1, b > 0, a < b. 

i=i 

In contrast to \2. 40\ ), the linear factors of w k do not depend k. 
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The first part of the following corollary gives an answer to Problem 3 in [4] , in the class 
of mean- field CFP's, while the second part recovers Proposition 1 in the above paper, 
in the aforementioned class of models. 

Corollary 2 For N large enough there do not exist mean-field Gibbs fragmentation 
processes on £l\m, with weights = (k — 1)!, k > 1 and = const, k > 1. 

Proof Recalling that = a^kl, k > 1, both assertions follow from fl2.50p which says 
that the asymptotics, as k — > oo of the two types of weights in question are not of the 
form required in Proposition 3. ■ 

Remark 6 In a recent paper FfOf a non mean-field Gibbs fragmentation process with 
weights Wk — (k — 1)!, k > 1 was constructed. The construction based on the Chi- 
nese restaurant model for simulation of uniform random permutation, results in a Gibbs 
fragmentation process with state transitions not obeying the mean field form 

• Spectral gap 

By virtue of (I2.43p . the spectral gap of the solvable CFP's considered is equal to the 
one of the Markov time homogeneous birth and death process \X^\t)\, t > with 
the rates of birth A r = A r ,jv = 0(1, 1){N — r), r = 1, . . . , N — 1 and rates of death 
jj r = ix r: N = ^-(^aN + roj, r = 1, . . . ,N. We shall employ Zeifman's method as 
described in [12], to find the spectral gap, say (3n, of the above birth and death process. 
Recalling that Xn = £*i = 0, consider the N — 1 quantities 

a r = a r ( 8 ) := \ r + /i r+1 - 8 r+ iX r+1 - -f, r = 1, . . . , N - 1, (2.63) 

where 5 — 5n — (5 r — 5 r ,N > 0, r — 2, . . . N — 1) is a vector of unknowns 5 r . The 
method states that 

(i) For any vector 5, 

min{a r , 1 < r < N — 1} < (3^ < max{a r , 1 < r < — 1} 

(ii) In the case of an ergodic birth and death process, there exists a unique vector 5, 
such that all N — 1 quantities a r are equal, so that their common value is equal to (3^. 
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In our case f!2.63j) conforms to 

a r = <p(l,l)(N-r) + ^(2aN + (r + l)b^-<j>(l,i)(N-r-l)S r +i- 
(r - 1) (2aN + rb) 

\- r = l,...,N-l. (2.64) 

Setting in (l2i)4j) S r = 1, r = 2, . . . , N - 1, we obtain 

a r = 0(1, 1) + aN + br, r = 1, . . . , N - 1, 

from which the following two -sided bound for the (3 is derived: 

0(1, 1) + aN + b < (3 N < 0(1, 1) + aN + b(N - 1). 

In particular, if b — 0, the preceding relation gives the exact value of the spectral gap 
(3 N = 0(l,l) + aiV. 
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